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ABSTRACT 


The  1988  Corn  Ear  Weight  Study  continued  to  analyze  the  forecast 
performance  of  models  estimated  using  surface  area  and  volume 
variables  to  predict  final  corn  ear  weight.  The  initial  study 
used  data  collected  in  Michigan  in  1986  and  1987.  The  1988  study 
was  expanded  to  Missouri  to  evaluate  the  performance  of  the  surface 
area  and  volume  models  in  a  corn  environment  that  is  more  drought 
prone  than  Michigan.  Two  models  based  on  diameter  measurements 
were  compared  to  models  estimated  using  the  operational  procedures 
from  the  Corn  Objective  Yield  Survey.  The  1988  Results  show  that 
the  research  models  have  mean  square  errors  that  are  40  to  57 
percent  lower  than  models  estimated  using  the  operational 
procedure.  These  results  support  the  findings  of  the  earlier  study 
which  indicate  that  the  research  models  have  a  superior  performance 
in  both  normal  and  drought  years.  Incorporation  of  a  surface  area 
or  volume  ear  weight  estimator  should  improve  the  new  Corn 
Objective  Yield  Models. 
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SUMMARY 


Forecasts  of  corn  yield,  acreage  and  production  are  obtained  from 
the  Corn  Objective  Yield  Survey,  which  is  conducted  by  the  National 
Agricultural  Statistics  Service.  Forecasts  of  ear  weight  and 
number  of  ears  are  used  to  predict  yield  per  acre.  The  purpose  of 
the  Corn  Ear  Weight  Study  was  to  develop  a  more  accurate  ear  weight 
estimator  using  surface  area  and/or  volume  variables  for  the  corn 
ear  in  addition  to  ear  length  measurements  made  in  the  operational 
program.  Diameter  measurements  made  one  inch  from  the  butt  of  the 
cob  and  two  inches  from  the  tip  of  the  cob,  and  length  measurements 
over  the  cob  were  used  to  compute  the  surface  area  and  volume 
variables . 

The  research  project  was  carried  out  in  Michigan  in  1986  and  1987. 
The  project  was  expanded  to  include  both  Michigan  and  Missouri  in 
1988.  Missouri  was  included  in  1988  to  determine  if  the  superior 
performance  of  the  research  models  in  Michigan  hold  in  a  corn 
environment  that  is  more  prone  to  drought.  1988  turned  out  to  be 
a  drought  year  in  both  states,  therefore  the  1988  results  provide 
valuable  information  about  the  ability  of  the  models  to  predict  ear 
weight  under  different  environmental  conditions. 

Two  research  models  were  estimated  using  surface  area  and/or 
volume  variables  in  the  1988  study.  These  models  were  also 
estimated  in  the  1986  and  1987  study.  Results  in  both  states  show 
that  the  research  models  have  mean  square  errors  that  are  40  to  57 
percent  lower  than  those  for  the  operational  models.  The  1988 
results,  which  support  findings  of  the  1986  and  1987  study, 
indicate  that  the  diameter  models  forecast  ear  weight  with  greater 
accuracy  in  both  normal  and  drought  years.  Therefore,  it  is 
recommended  that  the  new  diameter  measurements  be  added  to  the  Corn 
Objective  Yield  Program  in  the  1991  crop  year.  Data  will  need  to 
be  collected  for  a  few  years  in  order  to  build  the  models  before 
actual  implementation. 
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Forecasting  Corn  Ear  Weight  Using  Surface  Area 
and  Volume  Variables:  A  Final  Report 


Fatu  Wesleyl 


INTRODUCTION 


The  1988  Corn  Ear  Weight  Study,  a  continuation  of  the  1986  and  1987 
study,  examined  an  alternative  method  of  forecasting  corn  ear 
weight.  The  goal  of  the  study  was  to  develop  a  method  that 
produces  more  accurate  forecasts  of  ear  weight  than  the  operational 
method.  The  current  method  used  by  the  National  Agricultural 
Statistics  Service  (NASS)  forecasts  ear  weight  based  on  the  length 
of  the  ear  [4,5].  The  method  examined  in  the  Corn  Ear  Weight  Study 
forecasts  ear  weight  based  not  only  on  the  length  of  the  ear,  but 
more  importantly,  on  its  size.  Ear  size  is  estimated  by  computing 
the  surface  area  or  volume  of  the  ear.  The  1986  and  1987  results 
of  the  ear  weight  study  in  Michigan  showed  that  forecasts  derived 
from  the  research  models  were  closer  to  the  final  ear  weights  than 
forecasts  based  on  the  operational  model  [1]. 


The  project  continued  in  Michigan  in  1988  and  was  also  expanded  to 
Missouri  to  evaluate  the  performance  of  the  surface  area  and  volume 
models  in  a  corn  environment  that  is  more  drought  prone  than 
Michigan . 


Table  1  shows  NASS  official  corn  yield  estimates  for  Michigan, 
Missouri  and  the  U.S.  for  the  years  included  in  the  study.  The 
yields  for  1988  are  considerably  lower  than  the  1986  and  1987 
yields  because  1988  was  a  drought  year.  Therefore,  results  of  the 
1988  study  should  show  if  forecasts  based  on  size  are  still 
superior  to  forecasts  based  on  length,  during  a  drought  year. 
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Table  1.  NASS  official  corn  yield  estimates 
for  1986,  1987  and  1988. 


Year  State  Yield (bu/a) 


1986 

Michigan 

105.0 

Missouri 

116.0 

U.S. 

119.4 

1987 

Michigan 

97 . 0 

Missouri 

113 . 0 

U.S. 

119.8 

1988 

Michigan 

70.0 

Missouri 

76.0 

U.S. 

84 . 6 

The  Corn  Ear  Weight  Study  was  conducted  within  the  Corn  Objective 
Yield  Survey.  The  Corn  Objective  Yield  Survey  is  conducted  by 
NASS  to  forecast  corn  yield,  acreage  and  production.  This  monthly 
probability  survey  is  carried  out  from  August  to  November.  The 
survey  provides  forecasts  of  ear  weight  and  number  of  ears  which 
are  used  to  forecast  yield  per  acre.  The  Corn  Ear  Weight  Study 
examines  ear  weight,  one  of  the  components  used  in  the  yield 
forecast [6].  The  Corn  Objective  Yield  Manuals  [4,5]  and  the 
preliminary  report [1]  on  the  Corn  Ear  Weight  Study  give  a  detailed 
description  of  the  methods  used  to  obtain  an  ear  weight  forecast  in 
the  operational  program. 


In  this  report  "initial  study"  will  refer  to  the  1986  and  1987 
study  and  "final  study"  will  refer  to  the  1988  study. 
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METHODOLOGY 


Data  Collection 


Data  used  in  the  1988  study  and  in  the  initial  study  were  collected 
during  the  September  1  and  October  1  survey  periods  in  Michigan  and 
the  August  1  and  September  1  survey  periods  in  Missouri.  The  data 
were  collected  during  these  months  because  the  maturity  categories 
used  in  the  study  occur  during  these  months.  These  categories  were 
4  (milk) ,  5  (dough)  and  6  (dent) .  Maturity  category  3  was  excluded 
because  of  the  assumption  that  size  of  the  corn  ear,  at  this  stage 
of  maturity,  was  not  large  enough  to  show  a  significant 
relationship  to  the  final  ear  weight. 

In  the  initial  study  enumerators  examined  the  first  five  ears  of 
corn  next  to  a  designated  unit  to  determine  the  average  maturity  of 
the  corn.  They  measured  an  ear  in  the  unit  if  the  maturity 
classification  of  the  sample  field  was  1)  4,  5,  or  6,  based  on  the 
maturity  of  the  first  five  ears  of  corn  beyond  a  specified  unit  and 
2)  if  the  maturity  category  of  the  individual  ear  within  the  unit 
was  4,  5  or  6.  In  the  1988  study  an  ear  was  measured  if  the  sample 
unit  in  which  it  existed  met  the  first  condition,  which  is  the 
criterion  used  for  establishing  maturity  of  a  sample  unit  in  the 
operational  program.  It  is  assumed  that  eliminating  condition  (2) 
in  1988  did  not  create  inconsistency  in  the  data  collected  over  the 
years.  The  data  for  both  the  initial  and  final  studies  probably 
included  a  small  number  of  stage  3  ears  because  each  ear  in  the 
unit  was  not  examined  individually  during  the  final  study  and 
although  each  ear  was  examined  in  the  inital  study,  it  is  sometimes 
difficult  to  determine  the  maturity  of  an  ear  without  pulling  back 
the  husk. 

In  the  initial  study,  enumerators  made  three  measurements,  over 
the  husk,  on  ears  in  row  1  of  each  unit.  The  measurements  were: 
length  of  the  kernel  row,  diameter  of  the  cob  two  inches  from  the 
tip,  and  diameter  of  the  cob  one  inch  from  the  butt.  The  kernel 
row  length  measurement  was  made  using  a  tape  measure  and  the  butt 
and  tip  diameter  measurements  were  made  using  a  vernier  caliper. 
In  the  final  study,  the  kernel  row  measurement  was  replaced  with 
the  cob  length  measurement,  which  is  made  in  the  operational 
program.  The  replacement  was  made  because  the  kernel  row 
measurement  (over  the  husk)  did  not  show  a  more  significant 
relationship  to  ear  weight  than  the  cob  length  measurement,  which 
is  easier  to  make.  The  appendix  contains  questionnaires  used  for 
data  collection  in  1988. 
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Forecast  Model  Development 


In  the  operational  program,  the  ear  weight  models  are  estimated 
within  maturity  categories.  In  both  the  initial  and  final  studies, 
the  operational  and  research  models  were  estimated  across  maturity 
categories  because  of  lack  of  sufficient  data.  The  models 
estimated  are  described  below. 


Research  Models 


Research  models  for  forecasting  final  ear  weight  had  the  same  form 
in  both  the  initial  and  final  studies.  The  models  used  surface 
area,  volume  or  length  of  the  corn  ear  as  independent  variables. 
A  model  had  to  have  a  volume  and/or  surface  area  variable. 

Surface  area  and  volume  computations  used  a  diameter  measurement 
and  a  length  measurement.  The  diameter  measurements  were  made  one 
inch  from  the  butt  of  the  cob(BD)  and  two  inches  from  the  tip(TD). 
The  length  measurements  were  made  on  the  cob  over  the  husk(HL) . 

All  measurements  took  place  in  row  1  of  both  units,  except  for  HL 
which  is  an  operational  program  measurement.  HL  was  measured  in 
row  1  of  unit  2.  Sample  averages  were  computed  for  the 
measurements . 


The  definition  of  a  model  was 

Yi  =  +  B,X,i  +  B2X2i  +  ei 

where , 

Y^final  ear  weight  for  the  ith  sample  field, 

X^.  X2i  =  surface  area,  volume,  or  length 
variable,  and, 

e.  =  the  difference  between  final  ear  weight  and 
the  estimate  produced  by  the  model. 
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The  method  for  estimating  the  models  was  Least  Squares.  For 
inference  purposes,  this  method  assumes  that  Yi  has  a  normal 
distribution  with  mean  6o+B1X1+62X2  and  variance  a2  ;  and  that  e.  has 
a  normal  distribution  with  mean  0  and  variance  a2  . 

Two  research  models  were  estimated.  These  models  are  referred  to 
as  R1  and  R2 .  Both  models  were  estimated  in  the  initial  study. 
The  models  are: 


Model  Rl: 

Ear 

Weighti  = 

+ 

B.jSlj  +  B2BTi  +  e. 

Model  R2 : 

Ear 

Weighti  = 

*0  + 

6^2  f  +  B2BD.  +  e( 

Definitions  of  variables  included  in  the  models  are 


SI  =  HL*BD 

V2  =  (HL*TD2  )/4  *  7 r 

BT  =  BD*TD 


Where , 


BD  =  diameter  of  the  cob  measured  one  inch  from  the  butt, 

TD  =  diameter  of  the  cob  measured  two  inches  from  the  tip,  and 
HL  =  length  of  the  cob  measured  over  the  husk. 


Model  R1  includes  a  surface  area  variable,  SI,  and  an  interaction 
variable,  BT.  Model  R2  includes  a  volume  variable,  V2 ,  and  the 
butt  diameter,  BD. 


Operational  Model 


The  operational  model  is  a  weighted  average  of  two  regression 
models,  the  kernel  row  length (KL)  model  and  the  cob  length  (HL) 
model.  The  model  is  defined  as  follows: 


Let 

Y1i  =  “o+<MKL)i  +  < 

ekl  N  '  q2  kl) 

Y2,  =  B0+6,(HL),  +  ehli. 

ehl  f\J  N  (0 ,  a2  hl) 
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Then 


Y  ,•  =  ( a )  Y 1  j  +  (l-a)Y2,. 


(1) 


where  a=R2Y1i  /(R2Y1i  +  R2Y2l), 

and  R2Ylj  =  the  coefficient  of  determination 
for  the  kernel  row  length  model 

R2Y2i  =  the  coefficient  of  determination 
for  the  husk  length  model 


The  operational  model  referred  to  as  C_S  in  the  paper  on  the 
initial  study  [1]  is  referred  to  as  CURRENT  in  this  paper.  The 
model  estimated  in  (1)  differs  from  the  actual  operational  model  in 
that  the  actual  operational  model  is  estimated  within  maturity 
category  while  model  CURRENT  is  estimated  across  maturity 
categories  because  of  lack  of  sufficient  data. 


The  following  section  gives  the  methodology  used  to  compare  the 
research  and  operational  models.  The  research  models  are  not 
compared  to  the  new  operational  Corn  Objective  Yield  Models 
implemented  by  the  Statistical  Methods  Branch  in  1990.  This  issue 
is  addressed  in  the  recommendations  section  and  will  require  that 
the  new  diameter  measurements  be  collected  for  several  years  before 
they  can  be  incorporated  into  this  new  model  structure. 
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Model  Validation 


Model  validation  tests  were  conducted  on  the  research  and 
operational  models  in  order  to  compare  their  performance.  The 
tests  validated  the  accuracy  of  the  forecasts  and  the  stability  of 
the  model  coefficients.  Validation  for  forecasting  accuracy  was 
performed  because  the  primary  objective  in  developing  the  ear 
weight  models  was  to  obtain  more  accurate  forecasts.  The  stability 
of  the  model  coefficients  was  examined  to  determine  if  the  model 
relationships  are  consistent  over  years  and  therefore  reliable  for 
forecasting . 


Accuracy  in  Forecasting 


Two  methods  were  used  to  validate  the  models  for  forecasting 
accuracy.  The  methods  are  referred  to  as  Method  I  and  Method  II 
and  are  described  below.  Method  I  was  used  in  the  initial  study. 


Method  I 


Method  I  used  data  independent  of  the  estimation  data  for 
validation.  This  type  of  validation  gives  an  indication  of  the 
forecasting  ability  of  a  model  by  determining  how  the  model 
performs  using  data  other  than  that  used  for  its  estimation. 

The  validation  mean  square  error  (MSE)  was  used  to  compare  the 
forecast  performance  of  two  or  more  models.  This  mean  square  error 
is  computed  from  the  forecast  error  e.,  which  is  the  difference 
between  Y . ,  the  final  ear  weight  for  the  ith  sample  and  Yhatj ,  the 
ear  weight  forecast  for  the  same  sample.  Computation  of  the 
forecast  error  from  independent  (validation)  data  is  conducted 
using  coefficients  estimated  from  the  base  (estimation)  data  to 
predict  final  ear  weight.  MSE  can  be  expressed  as  follows: 

MSE  =  ( 1/n)  *  [  2  (e,- ) 2  ] 


Where , 

n  =  the  number  of  observations  in  the  validation 
data . 
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Figure  1  below  shows  the  data  used  for  estimation  and  validation. 
Missouri  had  data  for  only  one  year,  therefore,  a  random  sample  of 
1/3  of  the  data  was  used  for  validation. 


Figure  1.  METHOD  I:  Data  used  for  estimation  and  validation 
of  research  and  operational  models 


State 

Estimation  Data 
Peri od 

#  of  Obs. 

Val idation  Data 
Period 

#  of  Obs. 

MI 

1986  &  1987 

261 

1988 

113 

1988 

113 

1986  &  1987 

261 

MO 

2/3  of  1988 

64 

1/3  of  1988 

30 

Method  II 


The  second  method  of  validation  for  forecasting  accuracy  was 
recommended  by  NASS' s  Objective  Yield  Surveys  Section.  This  method 
compares  the  variances  of  the  models  but  makes  no  formal 
distinction  between  estimation  and  validation.  It  is  based  on  the 
premise  that  a  more  accurate  estimate  of  the  variance  of  a  model  is 
obtained  when  all  of  the  data  is  used  for  estimation  (none  is 
withheld  for  validation,  as  is  done  in  Method  I) . 


As  stated  earlier,  the  operational  model  is  a  weighted  average  of 
two  regression  models,  the  kernel  row  length (KL)  model  and  the  cob 
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length  (HL)  model.  The  following  procedure  is  used  to  estimate  the 
variance  for  this  model. 


Let 

Y1i 

=  <*0+<M  KL)f  +  ey1i. 

eYi  N  ( 0  >  °2  Y1 

Y2  j 

=  &0+fti  (HL)  ^  +  ey2j , 

eY2  A/  **  (  0  '  °2  Y2 

aY1,Y2  C°V  (eyi  f  ®y2) 

Y.=  (a)Yl,.  +  (l-a)Y21. 

Where  a=R2y1i  /(R2Y1i  +  R2Y2i)  / 

Then  Var(Y)  =[  (a2)a2y1  +  (l-a)2a2y2  +  2a(l-a)ay1  y2  ] 

ay1  Y2  is  estimated  by  Z.  rylirY2i/n 


Y1i= 

-  Y1  ■  and 

Y2i =  Y2i 

-  Y2,- 

Coefficient  Stability 


The  models  were  compared  for  cofficient  stability  by  using  an  F 
test  [3].  The  basis  of  the  F  test  is  a  comparison  of  a  model's 
error  sum  of  squares  when  the  coefficients  are  assumed  to  be  the 
same  to  the  error  sum  of  squares  when  the  coefficients  are  assumed 
to  be  different.  In  order  for  the  F  test  to  be  valid,  the 
variances  obtained  by  estimating  the  model  by  year  have  to  be 
equal . 


The  Barlett  test  for  equality  of  variances  [3]  was  used  to  test  the 
assumption  of  equal  variance.  The  appendix  describes  both  the  F 
test  and  the  Barlett  test. 
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RESULTS 


In  the  previous  section  the  methodology  used  to  estimate  the 
research  and  operational  models  were  given.  This  section  gives  the 
results  obtained  from  applying  the  methodology.  The  first  section 
gives  results  of  comparing  the  models  for  accuracy  in  forecasting 
and  the  second  section  gives  the  results  of  comparing  the  models 
for  coefficient  stability. 


Model  Performance 


Results  for  Methods  I  and  II  are  given  below.  In  Method  I,  a 
portion  of  the  data  is  used  for  estimating  the  models  and  the  rest 
of  the  data  is  used  for  validation.  Method  II  makes  no  distinction 
between  estimation  and  validation.  The  variances  of  the  models 
are  compared  when  all  of  the  data  is  used  for  estimation. 


Method  I 


Table  2  contains  estimation  statistics  for  the  research  and 
operational  models  which  were  computed  using  the  estimation  data 
given  in  Figure  1.  The  table  shows  that  estimation  mean  sguare 
errors  (mean  square  errors  from  the  regression  output)  of  the 
research  models  are  smaller  than  those  of  the  operational  model  and 
that  R2  s  of  the  research  models  are  higher  than  the  operational 
model's.  These  results  are  similar  to  results  found  in  the  initial 
study  [ 1  ]  . 


Table  3  gives  the  validation  MSEs  of  the  research  and  operational 
models.  The  table  shows  that  both  research  models  have 
substantially  smaller  mean  square  errors  than  the  operational 
models  for  all  validation  data.  The  reduction  in  mean  square  error 
when  the  research  models  are  used  is  at  least  40  percent.  These 
results  also  support  the  results  found  in  the  initial  study  [1]. 
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Table  2.  METHOO  I:  Estimation  statistics  for  Objective  Yield 
research  and  operational  models. 


State 

Year 

Independent  Variables 

Estimation 

MSE 

R 2 

MI 

1986  &  1987 

CURRENT 

KL 

.00503 

.2621 

HL 

.00473 

.3060 

R1 

(SI , BT ) 

.00222 

.6753 

R2 

(V2.BD) 

.00224 

.6730 

1988 

CURRENT 

KL 

.00554 

.3342 

HL 

.00407 

.5113 

R1 

(SI , BT ) 

.00255 

.6964 

R2 

(V2.BD) 

.00247 

.7058 

MO 

1988 

CURRENT 

KL 

.00950 

.2099 

ft 

HL 

.00792 

.3414 

R1 

(SI , BT ) 

.00440 

.6396 

R2 

(V2.BD) 

.00436 

.6430 

Table  3.  METHOO  I:  Validation  statistics  for  Objective  Yield 
research  and  operational  models. 


State 

Estimation 

Data 

Validation 

Data 

Model 

Variables 

MSE 

Percent  Reduction 
in  Forecast  Error 
R1  R2 

MI 

1986 

1988 

CURRENT 

(HL.KL) 

.0059716 

56.5 

57.0 

& 

R1 

(SI , BT ) 

.0026003 

1987 

R2 

(V2.BD) 

.0025675 

MI 

1988 

1986 

CURRENT 

(HL.KL) 

.0051602 

51.5 

48.7 

& 

R1 

(SI ,BT) 

.0025003 

1987 

R2 

(V2,BD) 

.0026473 

MO 

2/3  of 

1/3  of 

CURRENT 

(HL, KL) 

.0042620 

53.8 

51.4 

1988 

1988 

R1 

(S1.BT) 

.0019670 

R2 

(V2,BD) 

.0020713 
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Method  II 


Results  given  in  Table  4  are  for  Method  II,  which  uses  all  the 
data  to  compute  the  variances  of  the  models  and  makes  no 
distinction  between  estimation  and  validation.  The  table  shows 
that  the  research  models  have  substantially  lower  variances (43  to 
45  percent)  than  the  operational  model. 


Table  4.  METHOD  II:  Results  of  model  comparisons 


State 

Model 

Variables 

a1 

Percent  Reduction 

In  Variance 

R 1  R2 

MICHIGAN  2 

CURRENT 

(HL,KL) 

.00415 

44.6  45.1 

R 1 

(SI , BT ) 

.002299 

R2 

(V2,BD) 

.002284 

MISSOURI 

CURRENT 

(HL,KL) 

.00633 

43.4  43.3 

R1 

(SI , BT) 

.00358 

R2 

(V2.BD) 

.00359 

Michigan  analysis  was  conducted  using  data  for  three  years 
and  Missouri  analysis  was  conducted  using  data  for  one  year. 
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Coefficient  Stability 


This  section  gives  the  results  obtained  when  coefficients  of  the 
research  and  operational  models  were  compared  to  determine  if  they 
are  stable.  This  comparison  was  done  for  the  Michigan  models 
because  the  Missouri  data  was  for  one  year  only.  Table  5  shows  the 
coefficients  which  were  estimated  for  the  Michigan  models. 


Table  5.  Estimated  coefficients  for  the  research 
and  operational  models. 


COEFFICIENTS 


Model 

Variable(s) 

1986 

1987 

1988 

R1 

SI 

.000519 

.000386 

.000525 

BT 

.000104 

.000141 

.000134 

R2 

M2 

.000103 

.000115 

.000112 

BD 

.010303 

.009666 

.012002 

CHL 

HL 

.034228 

.044579 

.053529 

CKL 

KL 

.024564 

.032526 

.033862 

Determining  if  a  model's  coefficients  are  stable  involves  two 
steps.  First,  the  variances  of  the  model  over  time  have  to  be 
verified  to  be  statistically  equal  and  if  they  are,  an  F  test  is 
performed  to  determine  if  the  coefficients  are  stable. 

The  equality  of  the  variance  over  time  was  determined  by  conducting 
the  Barlett  test  [3],  which  is  described  in  the  appendix.  The 
Barlett  test  concludes  that  the  variances  are  equal  if  B,  the  test 
statistic,  is  less  than  \2 ,  the  (1-a)  percentile  of  the  Chi-square 
distribution  with  r-1  degrees  of  freedom.  Table  6  gives  the 
results  of  the  Barlett  test.  Values  of  \2  are  given  for  a=.05  and 
a=.10.  The  table  shows  that  B  is  less  than  J2  for  all  models 
except  the  operational  husk  length  model.  This  means  that  all 
models  except  the  husk  length  model  have  equal  variances,  therefore 
the  test  for  coefficient  stability  cannot  be  conducted  for  the  husk 
length  model. 
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Table  6.  Results  of  testing  research  and  operational 
models  for  equality  of  variance 


SSE 


Model 

Variable(s) 

1986 

1987 

1988 

B 

R1 

SI ,  BT 

.25258 

.30445 

.28040 

1.130 

R2 

V2,BD 

.26010 

.30380 

.27167 

.831 

CHL 

HL 

.48381 

.73263 

.45129 

8.856* 

CKl 

KL 

.53651 

.75625 

.61481 

3.880 

(2(a=.05)=5.99,  J2(ar=.  10)=4 

.61 

*  Statistically  significant 


Table  7  provides  the  results  of  the  F  test  for  coefficient 
stability  [3]  which  is  described  in  the  Appendix.  The  test 
concludes  that  the  coeffients  are  not  stable  if  F',  the  test 
statistic,  is  less  than  F,  the  (1-a)  percentile  of  the  F 
distribution  with  p  and  (nl+n2+n3-3p)  degrees  of  freedom.  The 
table  gives  F'  for  all  models  except  the  operational  husk  length 
model.  A  comparison  of  F'  to  F  shows  that  the  research  models 
exhibit  borderline  stability  while  the  operational  kernel  row 
length  model  is  clearly  unstable.  The  instability  of  the 
operational  model  coefficients  may  be  a  contributing  factor  to  the 
poorer  forecasting  ability  identified  previously. 


Table  7.  Results  of  testing  research  and  operational 
models  for  coefficient  stability 


SSE  SSE 


Model 

Vari able(s) 

Full 

Reduced 

F’ 

R 1 

SI  ,BT 

.83743 

.86485 

1 .99 

R2 

V2,BD 

.83557 

.86464 

2.12 

CHL 

HL 

1.66773 

1.7827 

- 

CKL 

KL 

1.9076 

2.0650 

5.02* 

F(a=. 

05)=2.60,  F(a=. 

.  10)=2. 13 

*  Statistically  significant 
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CONCLUSIONS 


The  1988  Corn  Ear  Weight  Study  continued  the  ear  weight  research 
which  occured  in  Michigan  in  1986  and  1987.  In  1988  the  study  was 
expanded  to  Missouri  to  evaluate  the  performance  of  the  surface 
area  and  volume  models  in  a  corn  environment  that  is  more  drought 
prone  than  Michigan;  however,  both  states  experienced  drought 
conditions  in  1988.  The  research  models  estimated  using  the 
Michigan  data  for  all  three  years  and  the  Missouri  data  for  1988 
have  mean  sguare  errors  that  are  4  0  to  57  percent  percent  lower 
than  those  for  the  operational  models.  Tests  performed  on  the 
Michigan  data  showed  that  the  research  models  have  coefficients 
that  exhibit  borderline  stability  while  coefficients  of  the 
operational  models  are  clearly  unstable.  These  results  support 
findings  from  the  1986  and  1987  Michigan  study  and  indicate  that 
the  research  models  predict  ear  weight  with  greater  accuracy  than 
the  operational  model  in  both  normal  and  drought  years. 


RECOMMENDATION 


Ear  diameter  measurements  should  be  added  to  the  Corn  Objective 
Yield  Program  beginning  with  the  1991  crop  year.  Research  indicates 
that  two  models  utilizing  these  measurements  perform  similarly  in 
predicting  the  final  ear  weight.  The  author  recommends  that  the 
surface  area  model  (Rl)  be  implemented  in  all  ten  corn  objective 
yield  states  because  it  has  a  slightly  simpler  form.  In  addition, 
research  should  be  conducted  to  determine  if  the  principle  of  using 
the  surface  area  of  the  ear,  instead  of  its  length,  improves  the 
performance  of  the  new  Corn  Objective  Yield  Models  [2]  implemented 
in  1990  by  the  Statistical  Methods  Branch.  With  this  new  approach, 
a  regional  model  sets  a  regional  estimate,  and  then  the  state 
models  are  constrained  to  produce  forecasts  that  are  consistent 
with  the  regional  model.  Data  will  be  needed  in  all  10  states  for 
at  least  5  years  in  order  to  build  the  regional  model.  Then  2 
additional  years  of  independent  data  should  be  used  for  model 
testing. 
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APPENDIX 


Stability  of  Coefficients 


Stability  of  estimated  coefficients  over  time  was  analyzed  by 
estimating  the  coefficients  by  year  and  using  a  statistical  test  to 
determine  if  coefficients  estimated  for  1986,  1987  and  1988 
Michigan  models  are  the  same  [3].  The  test  compares  the  error  sum 
of  squares  when  the  coefficients  are  assumed  to  be  the  same  with 
the  error  sum  of  squares  when  the  coefficients  are  assumed  to  be 
different.  The  larger  the  difference  between  the  two  sums,  the 
higher  the  probability  of  rejecting  the  hypothesis  that  the 
coefficients  are  the  same.  In  testing  the  equality  of  coefficients 
in  two  regression  equations,  each  of  which  includes  two  independent 
variables,  the  hypothesis  tested  was: 


Ho:  *01  *o2  *o3 '  *11  *12  *13'  *21  *22  *33 

Ha •  (*01^*02^*03)  and/ or  (Bn/B12/B13)  and/or  (B21/B22/B23) 


where , 


ol  • 

*11' 

B21  =  coefficients 
regression 

estimated 
equation , 

for 

and 

the 

first 

o2 ' 

*12' 

B22  =  coefficients 
regression 

estimated 

equation. 

for 

the 

second 

o3' 

*13' 

B23  =  coefficients 
regression 

estimated 
equation . 

for 

the 

third 

The  null  hypothesis  is  rejected  if  the  test  statistic 

[ SSE (R)  -  SSE(F)]/p  (Al) 

F '  = 


SSE (F) / (nl+n2+n3-3p) 


is  greater  than  F ( 1-a ;p, nl+n2+n3-3p) ,  the  1-a  percentile  of  the  F 
distribution  with  p  and  nl+n2+n3-3p  degrees  of  freedom.  The 
number  of  parameters  estimated  in  the  regression  function  is  p. 
The  probability  of  rejecting  the  null  hypothesis  given  the 
hypothesis  is  true  is  a.  In  (Al)  ,  SSE(R)  is  the  error  sum  of 
squares  obtained  by  assuming  the  coefficients  are  equal  (model 
called  reduced)  and  estimating  the  equation  all  three  years. 
SSE ( F)  is  the  error  sum  of  squares  obtained  by  assuming  the 
coefficients  are  different  (model  called  full) .  Computation  of 
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SSE(F)  involves  estimating  the  equations  for  each  year  separately 
and  adding  the  individual  sum  of  squares.  The  number  of  observa¬ 
tions  for  the  first,  second  and  third  years  are  nl,  n2 ,  and  n3 , 
respectively . 

The  test  for  determining  coefficient  stability  is  based  on  the 
assumption  that  variances  of  the  error  terms  from  the  regression 
for  the  first  year  (SSE.,),  second  year  (SSE2)  ,  and  third  year(SSE3) 
are  equal.  Since  there  are  more  than  two  variances  to  consider, 
the  Barlett  test  for  the  equality  of  variances  [3]  was  used, 
instead  of  the  F  test,  to  determine  if  the  models  satisfy  the 
assumption  of  equal  variance.  The  Barlett  test  concludes  that  the 
variances  are  equal  if 


B  <  X2  (1-a;  r-1) 


(A2) 


Where , 

B  =  (2 . 302585/C)  (nT-r)  (log10MSE  -  log1QGMSE) 

MSE=  the  arithmetic  mean  of  the  SSE's  weighted 

by  the  number  of  observations  for  each  year 

GMSE=  the  geometric  mean  of  the  SSE's  weighted 

by  the  number  of  observations  for  each  year 
nT=  the  sum  of  nl,n2,  and  n3 
r  =  the  number  of  years 

C  =  (1  +  (1/3  (r-1)  )  {  S;  (l/(n.-l))-(l/nT-r)  },  j=l  to  r 

/2(l-a;  r-1)  is  the  (1-a)  percentile  of  the  Chi-Square  distribution 
with  r-1  degrees  of  freedom. 
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UNITED  STATES  DEPARTMENT  OF  AGRICULTURE 
NATIONAL  AGRICULTURAL  STATISTICS  SERVICE 


PREHARVEST 


FORM  B-1:  CORN  YIELD  COUNTS  -1988 


Form  Approved 
OMB  Number 0535-0088 
Expiration  Date  7/31/89 
C  E  -1 20032B-1 


YEAR, CROP, FORM,  MONTH 

0-5) 

843 

| 

USE  FORM  B-1  FOR  ALL  PREHARVEST 
COUNTS  AND  MEASUREMENTS 


Date  ( _ 

Starting  Time  {military  time)  .  .  .  . 


Has  operator  applied  pesticides  with  organophosphorous  content  to  sample  field?  (  )  YES 
If  YES,  enter  latest  application  date _ and  name  of  pesticide _ 


370 


371 


(  )  NO 


UNIT  LOCATION  CODE 


1.  a.  First  visit  to  lay  out  unit . Code  1  « 

b.  Unit  relocated  this  month . Code  2  >  Enter 

c.  Sample  unit  laid  out  previously . Code  3  ^ 


Go  to  Item  3  when  coded  3;  otherwise  go  to  Item  2. 

ROW  SPACE  MEASUREMENTS 

2.  a.  Measure  distance  from  stalks  in  Row  1  to 

stalks  in  Row  2  . Feet  &  Tenths 

b.  Measure  distance  from  stalks  in  Row  1  to 
stalks  in  Row  5 .  Feet  &  Tenths 


UNIT  1 

UNIT  2 

302 

307 

UNIT  1 

UNIT  2 

303 

•  _ 

304 

•  _ 

305 

•  _ 

306 

STAGE  OF  MATURITY  Designated  Measurement  Areas: 


For  Month 

Use  Area  Bevond 

Maturity  Staqe 

Code 

Maturity  Staqe 

Code 

Aug  1 

Unit  1.  Row  1 

Pre-Blister.  .  .  . 

.  .2 

Dough  . 

.  .  5 

Sept  1 

Unit  1,  Row  2 

Blister . 

.  .  3 

Dent . 

.  .  6 

Octl 

Nov  1 

Unit  2,  Row  2 

Unit  2,  Row  1 

Milk . 

.  .  4 

Mature . 

.  .  7 

Husk  the  first  5  ears  or  silked  ear  shoots  beyond  the  unit  in 


the  designated  measurement  area  and  examine  for  maturity. 

If  ears  or  silked  ear  shoots  are  not  yet  present,  CHECK  (  ) 
and  complete  item  6  only. 

Ear  Number 

Total  of 

1 

2 

3 

4 

5 

5  ears 

3.  Maturity  stage  of  first  5  ears  or  silked  ear  shoots . 

301 

a.  Does  Item  3  have  3  or  more  Code  7  Ears? 


OR 

Will  this  field  be  harvested  within  the  next  3  days? 


▼ 


r  (  )  YES  -  to  either  question  -  Complete  B  -  2. 
(  )  NO  -  to  both  questions  -  Continue. 


TOTAL 

OF  5 
EARS: 

12  OR  LESS 

COMPLETE  ITEMS  6-9  ONLY. 

13-22 

COMPLETE  ALL  ITEMS. 

23  OR  MORE 

COMPLETE  ALL  ITEMS  EXCEPT  7  AND  8. 

NOTE:  Copy  maturity  of  ears  to  corresponding. cells  in  Item  4  on  back  of  form  (replace  any 
Code  2  ear  shoots  with  next  Code  3  ear  or  higher). 


FORM  B--1 :  CORN  (Continued) 


Ear  Number 


1 

2 

3 

4 

5 

4.  Enter  the  maturity  stage  of  each  of  the  first 

5  ears  Code  3  or  higher . 

320 

321 

322 

323 

324 

Enter  the  average  length  of  kernel  rows 

326 

327 

328 

329 

330 

(Item  4  ears) . 

.  Inches  &  Tenths 

• 

• _ 

•  _ 

•  _ 

• 

COUNTS  WITHIN  1 5  FOOT  UNITS 


6.  Number  of  stalks . 

7.  Number  of  stalks  with  ears  or  silked  ear  shoots 

(Item  7  cannot  exceed  Item  6  for  any  row) . 

8.  Number  of  ears  and  silked  ear  shoots 

(Item  8  MUST  equal  or  exceed  Item  7  for  any  row)  . 

9.  Number  of  ears  with  evidence  of  kernel  formation 

(Item  9  cannot  exceed  Item  8  for  any  row) . 


MEASUREMENTS  WITHIN  UNIT  2,  ROW  1 


UNIT  1  UNIT  2 


ROW  1 

ROW  2 

ROW  1 

ROW  2 

331 

332 

333 

334 

341 

342 

343 

344 

351 

352 

353 

354 

361 

362 

363 

364 

10.  Enter  the  number  of  ears  counted  in  Item  9  for  Unit  2,  Row  1  (Cell  363) _ . 

Measure  length  of  COB  for  all  EARS  counted.  The  number  measured  and  entered 
below  must  equal  the  number  counted  in  Cell  363. 


1. 

.  6. 

.  11- 

• 

2. 

•  7. 

.  12. 

• 

3. 

.  8. 

_ 13- 

• 

4. 

.  9. 

•  14. 

• 

5. 

•  10. 

_ 15. 

• 

16. 

.  21. 

•  . 

26. 

• 

17. 

.  22. 

,  • 

27. 

• 

18. 

23. 

• 

28. 

• 

19. 

1 

CM 

1 

• 

1 

• 

29. 

• 

20. 

25. 

• 

30. 

• 

Enumerator _ 

Did  your  supervisor  assist  you  in  working  this  sample? 
(  )  YES,  (  )  NO 


Total 

Length 

Total 

Ears 


Ending  Time 
( Military  Time) 

Enumerator 

Number 

Supervisor 

Number 


STATUS  CODE 


UNITED  STATES  DEPARTMENT  OF  AGRICULTURE 
NATIONAL  AGRICULTURAL  STATISTICS  SERVICE 

FORMB-R:  CORN  RESEARCH  STUDY -1988 


YEAR.CROP.FORM.MONTH 

(1-5) 


843 


Form  Approved 
O  MB  Number  0535-0088 
Expiration  Date  7/31/89 
C  E.-120032B-R 
Michigan 
Missouri 


1.  Does  Form  B  - 1,  item  3,  have  3  or  more  code  7  ears 

OR 

Will  the  field  be  harvested  within  3  days?  (  )  YES  --  Go  to  Form  B  -  2. 

(  )  NO  -  Continue. 


2.  Is  the  total  of  5  ears  given  on  Form  B  - 1,  Item  code  301, 18  or  greater? 

(  )  YES  -  Tag  the  3rd  and  4th  ears  in  row  1  of  each  unit  with  the  ID  tags  provided.  Complete  Item  3. 
(  )  NO  -  Return  to  Form  B  - 1 . 

MEASUREMENTS 


3.  Measure: 

The  length  of  the  cob.  Record  to  the  nearest  1/10  inch. 

The  diameter  of  the  ear  one  inch  from  the  butt  of  the  cob.  Record  to  nearest  1  millimeter. 
The  diameter  of  the  ear  two  inches  from  the  tip  of  the  cob.  Record  to  nearest  1  millimeter. 


UNIT  1,  ROW  1  UNIT  2.  ROW  1 


u 

Ear 

No. 

Cob 
Length 
(to  1/10  inch) 

Butt 

Diameter 
(to  1  mm) 

Tip 

Diameter 
(to  1  mm) 

U 

Ear 

No. 

Cob 
Length 
(to  1/10  inch) 

Butt 

Diameter 
(to  1  mm) 

Tip 

Diameter 
(to  1  mm) 

1 

1. 

2 

1. 

• _ 

.1 

2. 

2 

2. 

_ 

1 

3. 

2 

3. 

1 

4. 

2 

4. 

#  _ 

1 

5. 

2 

5. 

_ 

1 

6. 

2 

6. 

_ 

1 

7. 

2 

7. 

•  _____ 

1 

8. 

2 

8. 

•  _ 

1 

9. 

2 

9. 

•  _____ 

1 

10. 

• 

2 

10. 

•  _ 

If  more  than  1 0  ears,  continue  on  back 


FORM  B-R  CORN  (CONTINUED) 


UNIT  1,  ROW  1 


UNIT  2,  ROW  1 


u 

Ear 

No. 

Cob 
Length 
(to  1/10  inch) 

Butt 

Diameter 
(to  1  mm) 

Tip 

Diameter 
(to  1  mm) 

1 

11. 

•  _____ 

1 

12. 

1 

13. 

1 

14. 

m  _ 

1 

15. 

•  _ 

1 

16. 

1 

17. 

•  . 

1 

18. 

•  _ 

1 

19. 

•  _ 

1 

20. 

• 

1 

21. 

• - 

1 

22. 

1 

23. 

1 

24. 

#  _ 

1 

25. 

•  _____ 

1 

26. 

_ 

1 

27. 

•  _____ 

1 

28. 

•  — __ __ 

1 

29. 

•  . 

1 

30. 

•  _ 

U 

Ear 

No. 

Cob 
Length 
(to  1/10  inch) 

Butt 

Diameter 
(to  1  mm) 

Tip 

Diameter 
(to  1  mm) 

2 

11. 

2 

12. 

2 

13. 

2 

14. 

2 

15. 

2 

16. 

2 

17. 

•  

2 

18. 

2 

19. 

2 

20. 

• 

2 

21. 

2 

22. 

2 

23. 

2 

24. 

2 

25. 

2 

26. 

2 

27. 

2 

28. 

2 

29. 

•  

2 

30. 

•  _ 
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